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Abstract 

A Hamiltonian with two degrees of freedom is said to be super- 
integrable if it admits three functionally independent integrals of the 
motion. This property has been extensively studied in the case of two- 
dimensional spaces of constant (possibly zero) curvature when all the 

independent integrals are either quadratic or linear in the canonical 
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momenta. In this article the first steps are taken to solve the prob- 
lem of supcrintcgrability of this type on an arbitrary curved manifold 
in two dimensions. This is done by examining in detail one of the 
spaces of revolution found by G. Koenigs. We determine that there 
are essentially three distinct potentials which when added to the free 
Hamiltonian of this space have this type of superintegrability. Separa- 
tion of variables for the associated Hamilton- Jacobi and Schrodinger 
equations is discussed. The classical and quantum quadratic algebras 
associated with each of these potentials are determined. 
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1 Introduction 



A Hamiltonian system in classical mechanics with n degrees of freedom is de- 
scribed by a Hamiltonian function H{xi, . . . , Xn,Pi, ■ ■ ■ ,Pn) = H{x,p). The 
dynamics of such a system is described by Hamilton's equations 

OH . dH 

opi dxi 

The time rate of change of a classical observable i = i{x,p) is given by 

^ ^ " / a^a^ _ d^dH\ 

dt ' i^\dxidpi dpidxi) 
where { , } is the Poisson bracket. A Hamiltonian system is called "Liou- 
ville integrable" if it admits n functionally independent integrals of motion 
{Xi, . . . , Xn} which are mutually in involution, i.e. 

{X„X,} = 0, z,j = l,...,r2 (3) 

where one of these constants can be taken to be the Hamiltonian Q . The 
system is superintegrable if a further m integrals {Yi, . . . ,Ym,l < m < n—1} 
exist such that the set of constants {Xi = H, X2, . . . , X„, Yi, . . . , Y^} are 
functionally independent. The additional integrals have vanishing Poisson 
bracket with H, but not necessarily with each other or with the Xj's. A clas- 
sical Hamiltonian system is maximally superintegrable if m = — 1. There 
are then 2n — 1 functionally independent integrals of motion. The concepts 
of complete integrability and superintegrability have their analogue in quan- 
tum mechanics. In this superintegrable quantum mechanical system is 
described by n + m quantum observables {Xi = H, X2, ■ ■ ■ , Xn, Yi, . . . , Ym} 
which satisfy the commutation relations 

[H,X,] = HX,-X,H = 0, [H,Y,] = 0, [X„ 1^] = (4) 



where i,k = = l,..,m. For superintegrable classical Hamiltonian 

systems it is often the case that the elements of our set of constants are 
polynomial in the canonical momenta. The best known maximally super- 
integrable systems in Euclidean space En are the Kepler problem and the 
harmonic oscillator. All finite (bounded) trajectories in these two systems 
are closed. Moreover these are the only spherically symmetric potentials for 
which all finite trajectories are closed 

Systematic studies of superintegrable systems have been conducted for 
spaces of constant curvature in two and three dimensions ^ |^ ||, ^ |T0[ . 
In particular, a complete classification all superintegrable systems in the real 
Euclidean spaces E2 and £"3 with at most second order integrals of motion was 
given ^ H . More recently a relation between superintegrable systems 

and generalized Lie symmetries has been established [0, as well as their 
relation to exactly solvable problems in quantum mechanics 0]. Recently 



13, 14, Q5[ it has been possible to classify all maximally superintegrable 



systems for spaces of constant curvature (possibly zero) in two dimensions 
for which all the extra constants of the motion are at most quadratic in the 
canonical momenta. 

A natural question to ask is whether the concept of superintegrability 
is restricted to spaces of constant curvature. The purpose of this article is 
to show that this is not so and to start a study of superintegrable systems 
in more general Riemannian, pseudo-Riemannian and complex Riemannian 
spaces. More specifically, we consider real two-dimensional spaces and search 
for Hamiltonian systems allowing additional constants of the motion that are 
at most quadratic in the momenta. 
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To make initial progress on this problem we first need to know which 
Riemannian spaces in two dimensions have associated with them more than 
one classical quadratic constant of the motion. This is a problem that has 



been comprehensively solved by Koenigs in a note written in volume IV 
of the treatise of Darboux . 



In addition to being of intrinsic interest, additional motivation for this 
problem comes from the observation that all two-dimensional Riemannian 
spaces can be embedded in the three-dimensional Euclidean or pseudo- Euclidean 
space. Consequently any such two-dimensional classical motion is equivalent 
to a constrained motion in three dimensions. It is also possible to interpret 
the motion, via general relativity, as motion in a two-dimensional gravita- 
tional field. 

Given that we have a Riemannian space in two dimensions with infinites- 
imal distance 

ds^ = gij{u)du^du\ i,j = 1,2 (5) 
and u = {u^, u^), the classical Hamiltonian has the form 

H = ^g''p^PJ + V{u) (6) 

and the corresponding Schrodinger equation can be taken to have the form 

= -^9^' {y/^g^'d^.^) + V{u)^ = (7) 

where g = det{gij). For the classical Hamiltonian H our problem is to look 
for potentials V{u) and Riemannian spaces specified by the metric gij for 
which there are at least two extra functionally independent constants of the 
motion of the form 

Ai = a''iu)pipj + h{u) (8) 



or 

A2 = a%u)pi + c{u) , (9) 

other than H. One well known way of solving the corresponding classical 

problem is to use Hamilton-Jacobi theory. The crucial equation to solve is 
then the Hamilton-Jacobi equation obtained from the equation H — E via 
the substitution Pi — ■§^, that is 

This equation is sometimes solvable by the method of separation of variables 
using the additive separation Ansatz 

S ^ Si(u\a,E) + S2{u^,a,E) . (11) 

The corresponding Schrodinger equation can also be solved by separation of 
variables with the product Ansatz 

^^Mu\\E)i;2{u^X,E). (12) 

The quantities Aj are constants of the motion if 

{Xi,H}^0. (13) 

For A2 this implies that a'^{u) is a Killing vector and a\u)pi is a symmetry of 
the free Hamiltonian [H without V{u)). In the case of Ai this imphes that 
a^^{u) is a Killing tensor. Such tensors are directly related to the notion of 
additive separation as described above. We note that for constants of the 
type A2, the condition implies c{u) — 0. It is also clear that for every constant 
linear in the momenta, its square is a constant quadratic in the momenta, 
that is, of the form of Ai. 
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As mentioned earlier, Darboux and Koenigs have given a comprehensive 
analysis of when a two-dimensional Riemannian space admits more than one 
quadratic constant. In Section ^we summarise some of these results [|I8| , |19 



In the remaining sections we concentrate on a particular space with a Killing 
vector and two Killing tensors. Section ^ deals with the free Hamilton- Jacobi 
equation and we show that the Schrodinger equation allows separation of 
variables in three different coordinates systems which we determine explic- 
itly. Potentials that allow separation of variables in these systems are then 
introduced. In Section ^ we find all potentials with this superintegrability 
property. We then discuss in Section |^ the various surfaces that may be 
represented by the infinitesimal distances that we have and the consequent 
special functions that arise from the corresponding Schrodinger equation. 

2 On geodesies with quadratic integrals 



In 1889 G. Koenigs [0 wrote a note in the last volume of Darboux's treatise 
"Theorie generale de surfaces" , the title of which coincides with the title of 
this section. This note contains a summary of results which are the solution 
of the problem outlined in the Introduction, viz. when does the free Hamilto- 
nian of a two-dimensional Riemannian space admit more than one quadratic 
constant of the motion. The analysis was performed over the field of com- 
plex numbers and must be modified over the reals. What Koenigs did was 
to write the infinitesimal distance for a general two-dimensional Riemannian 
space in the form 

ds^ = Af{x^y)dxdy . (14) 
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This can always be done in two dimensions over C. The corresponding free 
Hamiltonian then has the form 

^= 27(b) 

By making the requirement that there is a second order Kilhng tensor of the 
form 

X ^ a'^ {u)piPj , (16) 
Darboux and Koenigs estabhsh the foUowing propositions. 

1. Any two-dimensional Riemannian space that admits more than one 
Kilhng vector must be a space of constant curvature and admit three 
linearly independent Killing vectors. 

2. Any two-dimensional Riemannian space that admits more than three 
Killing tensors is a space of constant curvature. It then actually ad- 
mits five linearly independent Killing tensors which are all bilinear ex- 
pressions in the Killing vectors. The sixth bilinear combination is the 
Hamiltonian itself. 

3. Any two-dimensional Riemannian space that admits precisely three 
linearly independent Killing tensors will be a Riemannian space of rev- 
olution. In fact there will be one Kilhng vector and two Killing tensors. 

Two-dimensional Riemannian spaces of this latter type were distinguished 
to be of four types. The infinitesimal distances of these types are given by 

(I) ds^ — {x + y)dxdy. 



(III) ds^ = (ae-^ + 6e-^-^) dxdy. 

/ x-y , 

a e 2 +62 +6 

(IV) ds^ = — ^ dxdy. 

I e 2 — e 2 j 

It is the first of these infinitesimal distances that we analyse in some detail 
in the next section. We shall call the spaces "Darboux spaces" and denote 
then by Di, D2, D3 and D4 respectively. 



3 The free particle and separation of vari- 
ables in a Darboux space of type one 

If we consider the first space of Darboux's list and look at real forms of this 
space only, it is convenient to make the new choice of variables 

X — u + iv, y — u — iv . (17) 

The corresponding infinitesimal distance can then be taken as 

ds^ = 2u{du^ + dv^) , (18) 

and the corresponding Hamiltonian has the form 

H=l;ipl + Pl)- (19) 

Associated with this Hamiltonian are three integrals of the free motion, two 
quadratic and one linear. 

K ^ 

Xi = PuPv-^{pI+pI) 

X2 = Pv{vpu-up^) - —{pI+pD . (20) 
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These three integrals satisfy the polynomial Poisson algebra relations, 



{K,X,}^2H, {K,X2}^-Xi, {X,,X2}^2K^. (21) 
They cannot be functionally independent and in fact satisfy the relation 

+ Xf + K'^^0. (22) 
For the analogous quantum problem it is sufficient to consider the operators 

K — —idv 

X2 = -l[dv,vd^-ud,]^ + ^{dl + d',) (23) 

where [^4, B]+ = AB + BA. The quantum versions of the quadratic constants 
are obtained via the formula 

A = --^a, {a^'Vgdj) . (24) 

These operators have the same commutation relations as for the classical 
constants with the Poisson bracket replaced by the commutator bracket. 

[k, Xi] = 2iH , [k, X2] = -iX^ , [Xi, X2] = -2ik^ . (25) 

There is also the operator relation 

4HX2 + X^ + k^^O. (26) 

The question we address in this section relates to the various possible ways 
that separation of variables can be achieved in the case of free classical motion 
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or its quantum analogue, the free Schrodinger equation. The criteria for this 

to occur is the same in either case. Classically, if we have a general quadratic 
first integral A and free Hamiltonian 

H = ]^g'^{u)p,p,, (27) 

and if the characteristic equation, 

0, (28) 
has two distinct roots pi and p2, the Hamiltonian will have Liouville form 



H = y^^^>fp2 (29) 



when written in terms of the new variables pi, P2- That is 

^(pi) Ppi + r{p 2)pl 

Pi + P2 

In this form, both classical and quantum systems can be solved by the sep- 
aration of variables Ansatz. 

If we want to classify all different separable coordinate systems for a given 
Hamiltonian we need to know how many essentially different quadratic first 
integrals are possible. To decide on the notion of equivalence we first observe 
that the variable v does not explicitly appear in the metric tensor, that is, 
it is an ignorable variable. This means that the transformations v ^ v + b 
form a one- dimensional Lie group. Accordingly, we determine the notion of 
equivalence to mean that two quadratic integrals are equivalent if they are 
related by a motion of this group. Consequently the most general quadratic 
constant can be written 

A = aXi + bX2 + cK^ (30) 

to within the addition of a multiple of H. The second order elements Xi 
transform under the adjoint action according to 

Xi — > exp{aK)Xiexp{—aK) 
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e:x.p{aAd{K))Xi 



X, + a{K, Xi] + ]^a\K, {K, X,}} + • • • (31) 



or specifically 



Xi ^ Xi + 2aH 

X2 X2-aXi-a^H . (32) 

There are three classes of possible quadratic first integrals under this equiv- 
alence relation. Typical representatives are 

Xl + a^^^ X2 + ai^^ (33) 

We can now explicitly demonstrate the separable coordinates in each of these 
cases. 

1. Separating coordinates associated with Xi + aK^ 
If we choose a representative to be 

L = J5^l + sinhcX^ (34) 

the corresponding roots of the characteristic equation and hence new 
variables are 

r = pi = -2{Cu + v) 

s^p2 = ^{u-Cv), C^e-". (35) 
In terms of these coordinates the Hamiltonian has the form 

(36) 



C(s-r) VC2 

and the corresponding quadratic constant in these coordinates is 



^ = ^-C{^{.C-^ri + sri]. (37) 
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2. Separating coordinates associated with X2 + aK^ 

Taking the second representative in the hst (^), that is L = X2 + aK'^, 
a convenient choice of new variables ^, 77 are related to the roots pi by 

p, = 7^\2a-r]^), p2 = -e'(2a + a. (38) 

The corresponding classical Hamiltonian then has the form 

2(e + ^')(e-^' + 2a) ■ ^ ' 

The associated constant of the motion in the new coordinates ^ and rj 
is 

^ rj\2a-v^)pl-e{e + 2a)pl 

The defining coordinates u, v are written in terms of the new coordi- 
nates ^, rj via 

u = ^ie-v') + a, v = ir] (41) 
which looks like displaced parabolic coordinates in the u, v plane. 

3. Separating coordinates associated with 

For the last representative, , we need only the coordinates u^v and 
to recognise the fact that K = p^. 

We conclude this section by discussing the solutions to the free particle 
and free Schrodinger equation in these three cases. 

In case 1 above it is more convenient to choose the variables according to 

u = r cos 9 + s sin 6* , v = —r sin 9 + s cos 9 . (42) 

The classical Hamilton- Jacobi equation then has the form 

H = y^'^ = E (43) 

4(rcos6' + ssm6') 
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which has the general separable solution 

6£/ cos t) oE sm fc' 

The corresponding free Schrodinger equation 

= -77 ^(^r + 5')* = (45) 

4(r cos 6* + s sm 6*) ^ r- v / 



has the typical product solutions 



3/2N 



AEcosOJ V 4£;sin^y ^ ^3 V 4£;cos^ 

X Ci ( -\/4Esin^ fs + ^y^^ 1 , (46) 

5 V 4Esin^; y ' ^ ^ 

where C,^(2;) is a solution of Bessel's equation. 

In the second case the classical Hamilton- Jacobi equation is 



TT _ \di J \dri / _ rp /^-r^, 

and has a general solution of the form 

J ^J2E^^^2E^^^^ di + j ^-2Er}^ + 2Ecr/2 + A d-q , (48) 

which can be expressed in terms of elliptic integrals. The corresponding 

Schrodinger equation has a solution of the form \1' = ipi{i)ip2{'n) where the 
il^i satisfy 

(a| + 2^e' + 4^cf + A)V^i(0 = 

{d^^ - 2Er]^ + AEcT]'' - X)Mv) = 0. (49) 

These equations are readily identified as the equations for the anharmonic 
oscillator. 
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In the third case the classical Hamilton- Jacobi equation is 

which has separable solutions 

S^l-{^Eu-k''f''' + kv. (51) 
The separable solutions to the corresponding free Schrodinger equation 

-l^i^u + dD^^E^ (52) 



have the form 



3 \ 3 I AE 



'^ = \h-irBCr\-^AE\u- — \ |e-V (53) 



It is clear that the actual solutions to the classical motion or the corre- 
sponding Schrodinger equation depend on the range of values assumed by 
the various real variables, that is, on exactly which real manifold we are 
considering. 

4 Integrable and superintegrability systems 
for the Darboux space of type one 

In this section we address the problem of superintegrability for the Hamilto- 
nian 

H^^{pI + pI), (54) 
that is, look for potentials V{u, v) for which 

H ^ H + V{u,v) (55) 
15 



admits at least two extra quadratic integrals. The way to solve this problem 
is as follows. First we consider that we already have one quadratic first 
integral 

L = a{u, v)pI + h{u, v)puPv + c{u, v)pl + d{u, v) . (56) 

We know that the quadratic part of L (i.e. that part obtained by putting 
d{u, v) = in (p6D) must correspond to one of the three possibilities outlined 
in the previous section. For each of these possibilities separation of variables 
is possible in coordinates a, [3 where u = u{a,P), v = v{a,(3). The addition 
of a potential implies that separation is preserved. As a consequence of this 
H can be written as 

a{a) + T{fj) 

and the corresponding first integral will be 

a{a){pl + g{f3))-T{P){pl + fia)) 

a(a) + r(/5) • ^^^^ 

The next step is to impose the condition that there is a further quadratic 
first integral and see what conditions this imposes on the functions /(a) and 
g{l3). If we do these calculations systematically we arrive at the following 
three cases. 



1. 



pI+pI , bi{4u^ + v^) ^2 , b. 



^ ^ ru_j:^ ^ + - + ^ . (59) 

Au Au u uv^ 

The additional constants of the motion have the form 



Rl = X; 



2 



4tt u v^u 



i?2 = K^ + hiv^ + ^ (60) 
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and the corresponding quadratic algebra |T6|, relations are deter- 
mined by 

{R,Ri} = 8HRi + 6Rl + 16b2R2-32bibs 
{R,R2} = -8HR2-l6biRi 

R^ = -16HRiR2-4:Rl-16b2Rl-6^b3H^ -IGbiRl 

+646i63i?2 + 2566162&3 (61) 

where R = {Ri, i?2}- The Hamiltonian clearly separates in the coordi- 
nates u and V as well as the coordinates ^, r] given by w = jiC'^ — i]'^) +a, 



V = ^rj. This can be seen from the explicit form 



H 



+ 



462 + ^ 



2(^2 _ ^2 + 2a) 



(62) 



The corresponding quadratic quantum algebra relations are 



[R,Ri] 
[R,R2] 
R' 



-6Rl - 8HR1 + 1662^2 + 26i(3 + I663) 
8HR2 - IGbiRi 



+AR^. - 8H 



Ri,R2 - m2Rl - IQbiRl 



(63) 



-46i(ll + 1663)^ - 4(3 + 1663)//^ + 166162(3 + I663 



where R = [Ri, R2] 



2. 



H 



pI+pI ^ ^ ^ ^ ^ Q3(^^ + v'^) 
iu u u u 



The additional constants of the motion have the form 



Ri 
R2 



2aiv 2a2{u^ — v"^) 2a3v{u'^ — v"^) 



+ 



u u 



u 



(64) 



(65) 
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and the corresponding quadratic algebra relations are determined by 
{R,Ri} = -8H^ + 16a3R2 + 8{al + iaias) 
{R,R2} = l6a2H-16a3Ri 

R^ = 16H^R2-16a3Rl + 32a2HRi-16a3Rl 

-16{al + 40103) i?2 - 64aia2 (66) 

If we change the coordinates according to m = r cos 6* + s sin ^, v — 
—r sin 9 -\- s cos 9 the Hamiltonian assumes the form 

^ _ + + 4ai + 4a2(-r sin 6' + s cos 9) + 4a3(r^ + s^) 

4(r cos 6' + s sin 9) 

which clearly also separates in these coordinates. 

The commutation relations of the corresponding quantum algebra are 

[R,Ri] - 16a3i?2 + - 8(0^ + 4aia3) 
[^,^2] = -1603^1 + 1602^ 

R^ = -1603^^-1603^^ + 16^2^2 + 3202^1 

-16(o^ + 40103)^2 + 64(o^ - oio^) . (68) 

3. The third potential gives rise to a Hamiltonian of the form 

iy = ^ + ^. (69) 

There are three extra constants associated with this Hamilonian, 

R^^Xi- — , i?2-^2-— and K. (70) 

u u 

The associated Poisson bracket relations are 

{ir,i?i} = 2H 
{K.R^} = -i?i 

{RuR2} = 2K{K^ + 2a) (71) 
18 



and the corresponding functional relation amongst these constants is 

AHR2 + Rl + K'^ + AaK^ = . (72) 

The commutation relations associated with the corresponding quantum 
problem have the form 

[K,Ri] = 2iH 
[K,R2\ = -iRi 

[RuRi] = -2ik{k^-2a) (73) 
and the identity amongst the defining operators is 

AHR2 + Rl + k*- 4ak^ = . (74) 

Upon examination of the various superintegrable potentials we have con- 
structed we see that by multiplying the equation H — E hy a suitable factor 
we essentially recover a variant of one of the superintegrable systems already 
classified for spaces of constant (or zero) curvature. For the first potential 
above, the equation H — E may be written 

pI + pI + ^i(4w' + v^) + 462 + ^ - 4£;i^ = . (75) 

V 

This equation is known to have separable solutions in coordinates it, v and 
associated parabolic coordinates ^, 77 given hy u — \{^^ — rf), v — ^rj. With 
the second potential, H — E becomes 

pI + pI + 4a3 {u^ + 1;^) + 4ai + 4a2i; - 4:Eu = (76) 

and the third, 

pl + pl-4:Eu + 4:a^0. (77) 
19 



This observation is crucial to the whole programme that we will undertake 
which aims at finding all superintegrable systems associated with a curved 
space in two dimensions and having quadratic constants. 

All three of the above systems are special cases of the superintegrable 



systems found in E2 @, |TT|. They were shown to be exactly solvable in Ref. 
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5 Embeddings of a Darboux space of revolu- 
tion of type one 

It is clear that the infinitesimal distance 

ds'^ = 2u{du^ + dv^) (78) 

does not uniquely determine a manifold. This then gives rise to the question 
of just what sort of surfaces can this infinitesimal distance represent. A 
particular choice of such a surface would determine the range of variation 
of the parameters u, v which in turn enables the solution of the geodesic 
equations in the case of classical mechanics and the quantum mechanics of a 
point particle. It is known that any two-dimensional Riemannian space can 
be embedded in a three-dimensional Euclidean space of indefinite or definite 
signature. In this section we look at a number of natural embeddings and 
discuss their associated geodesies and quantum mechanics. The infinitesimal 
distance that we are dealing with can be embedded in three-dimensional 
Euclidean space via the formulas 

X = V2ucosv, Y = V2usmv , (79) 
20 



where u > ^, Vq < v < 27r + Vq, sin f — ^/2u + 1 and F{(f, k) is an elliptic 
integral of the first kind. This embedding gives the infinitesimal distance 

dX^ + dY^ + dZ^ = 2u{du^ + dv^) . (81) 

To do quantum mechanics on this surface let us first look for separable solu- 
tions to the free Schrodinger equation. A typical solution has already been 
found in the previous section, viz. 

where m is an integer. As it > | and we see that u — ^ is not a singular 
point of the separable equation in u, we can impose a condition of the form 

a^Q,^;)+6^„Q,^;)=0 (83) 

together with the periodic boundary condition \£'(w, v) = ^(?i, v + 27i), which 
is already satisfied. If we take o = 1, = then E > 0, otherwise there is 
no solution satisfying the boundary condition at u — ^. U E > then we 
can find a suitably behaved solution that vanishes as u — > oo and satisfies 
the boundary condition at m = |, viz. 

* = {UU'f' {J.{U)J_.{U') - J.{U')J-l{U)) (84) 

where U = |a/4E (u - ^f^^ and U' = IVIE (| - These solutions 

are the analogue of the scattering states on this manifold subject to the 
boundary condition we have adopted. 
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An interesting embedding in pseudo-Euclidean space is given by 
X = V2uv 




for which dX'^ + dY^ - dT'^ = 2u{du^ + dv'^). In this case the variables vary 
over the ranges — cxD<t)<oo,0<M<cxD. We could indeed do an analysis of 
the free Schrodinger equation on this surface and come to a similar conclusion 
if we imposed the condition that the wave function is zero at u = 0. However 
if we consider the first potential (|59D and choose bi = — h = \ {\ — 7^) 
for real (3 and 7 > and if we write the solutions to Schrodinger's equation 
in the form \E' = U{u)V{v) then two independent solutions of the separation 
equation satisfied by V can be taken as 

If we wish to interpret these solutions as being associated with an angle 
variable which varies in the range < vq < v < Vq + 2tt and then we would 
require the periodic boundary conditions 

Vivo) = Vivo + 27T) 

V'ivo) = V'{vo + 2n). (87) 

The possibility of imposing these boundary conditions depends on whether 
V = occurs inside the domain of v. If it does not then the spectrum is 
determined from the condition 

W [V^ix) -V^ix + 2n), V4x) - V_(x + 2n)] |,^,^ = . (88) 
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If 1) = is included then the same conditions no longer work as this is a 
regular singular point of the equation. Indeed if t^o = and we assume 
7 > |, then we choose the solution V+ and impose the condition 

V+{27r) = (89) 

as V^+(0) is already zero. The quantisation condition is then determined by 

iFi(^^(l + 7)-^,l + 7,4/37r2^ =0. (90) 

For the u separation equation the range of variation of the variable u > | is 
clear and u = | is not a singularity of the the separation equation. We can 
accordingly take typical solutions to be 



where = + 4^2 ^ A*) ~ \- To obtain a solution that vanishes as 

u oo requires that 02 = 0. The remaining boundary condition becomes 



(92) 



This condition determines the nature of the discrete spectrum. For large 
eigenvalues the discrete spectrum is given by 



E ^ -2^(3^n (93) 

for suitable large integer n. 

If we consider the second potential ( |64D then putting 03 = —a^ the equa- 
tion for V{v) has solutions of the form 

V = diDJ2,/^(v-^]) +d2DJ-2,/^(v 



diV+ + d2V_ (94) 
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where 

1 



4/x+^)-^ (95) 



As there are no singularities in v in this equation and we can require that 

Wq < V < Wq + 2% 

V{wo) = Viwo + 27i) 

V\wo) = V'{wo + 27r) (96) 

which is equivalent to 

W [V+{x) - V+{x + 27r),V4x) -V4x + 2^)] = . (97) 
The solutions for the function U (u) that are well behaved for large u are 

U{u) = D,(2V^(u-^^^ (98) 
where p = ^ (4ai - 4/i + J) - |. 

6 Conclusion 

In this article we have examined one of the four spaces of revolution listed by 



Koenigs |T^. For the space that we have considered, it has been shown that 
there are three potentials that can be added to give superintegrable Hamilto- 
nian systems of the type we seek. In each of these cases we have exhibited the 
various inequivalent ways in which a separation of variables can be achieved 
for both the classical and quantum equations that result. This is equivalent 
to determining the various inequivalent ways in which a Hamiltonian can 
be written in Liouville form ( p7D for suitable separable coordinates a, (]. In 
particular we note that each of the three superintegrable systems we have 
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examined are such that when we write out the classical equation H = E 
and factor out the denominator we recover a variant of a superintegrable 
system corresponding to flat space [H]. This is an example of what is called 



coupling-constant metamorphosis |2^. It has been proven in [T^ that all 
of the superintegrable systems in the plane are such that the bound states 
energies can be calculated algebraically. In all cases the Hamiltonian lies in 
the enveloping algebra of s/(3,]R). We conclude that analogous statements 
apply to the superintegrable systems that we have found. 
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